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Abstract 

This paper presents an analytical approach of the propagation of an 
acoustic wave through a normally distributed disordered lattice made up 
of Helmholtz resonators connected to a cylindrical duct. This approach 
allows to determine analytically the exact transmission coefficient of a 
weakly disordered lattice. Analytical results are compared to a well- 
known numerical method based on a matrix product. Furthermore, this 
approach gives an analytical expression of the localization length apart 
from the Bragg stopband which depends only on the standard deviation 
of the normal distribution disorder. This expression permits to study on 
one hand the localization length as a function of both disorder strength 
and frequency, and on the other hand, the propagation characteristics on 
the edges of two sorts of stopbands (Bragg and Helmholtz stopbands). 
Lastly, the value of the localization length inside the Helmholtz stopband 
is compared to the localization length in the Bragg stopband. 

Wave propagation in random media, waveguide, scattering of acoustic waves. 
PACS ll.80.La ; 42.25.Dd ; 43.20.Mv ; 43.20.-f 

1 Introduction 

Since 1958 and the work of P. W. Anderson on the locaHzation of wave in a 
random media [T], the propagation of waves in complex media has been at the 
center of many works. Number of them have been undertaken on the effect of 
disorder on the wave propagation and many research fields have been concerned. 
The first studies appeared in sohd state physics [SJ |3j and in the propagation 
of EM waves in random medium [U [5] . The propagation of classical waves in 
disordered media [51 El [5] and the localization of elastic waves [HI UHl E] have 
been also dealt with theoretical and experimental studies and the applications, 
for example in acoustics [HI [T31 HH [H] or in geophysics [HI [T7] are numerous. 
In the same time, number of works have proposed experimental results in ID 
[nHHdHldl], 2D [HHD] and 3D media [12 Hg. 

Among these numerous references, only a few analytical studies of Ander- 
son localization have been published. First of all, Kholer et al studied, in the 
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70's, the wave propagation in a one-dimensional medium with random index of 
refraction or a transmission hne with random capacitance per unit length using 
the radiative transport theory In this study, the "disorder" is measured 

by a time-invariant refractive index field. The frequency dependence of the lo- 
calization length for acoustic and electromagnetic waves in a one-dimensional 
randomly layered media is also studied analytically [53] . An analytical theory of 
a pulse propagating in a one-dimensional layered media [25] is also established 
by using asymptotic methods for stochastic differential equations [IS] and a 
low frequency limit is considered to study the localization of elastic waves in a 
plane-stratified media [57] . 

Secondly, the problem of pulse backscattering from a randomly stratified 
media is considered where the wave speed fluctuations depend on time and on the 
range coordinate. The disorder is measured by a time-varying refractive index 
field and the time variation is parametric [28j . Radiation transport equations 
are used to describe the propagation and the wave localization in time domain 

m- 

In the same way, some works present analytical results by using approximate 
methods like the Coherent Potential Approximation [5] [6] , in the case of small 
amount of impurities [U [3D] , for the low frequency case [3T] or for asymptotic 
behavior [55] . 

Nevertheless, the Transfer Matrix method [ni[I3] is generally used to simulate 
the wave propagation in a random lattice |33| . 

The present paper proposes a new analytical approach to study the prop- 
agation of acoustic waves in a ID random media. This method is based on 
a recursive relation describing the wave propagation. It uses the properties 
of the normal distribution disorder to propose an analytical expression for the 
transmission coefficient depending on the standard deviation of the disorder 
distribution and on the number of lattice cells. 

The section II presents a general study of the wave propagation in a ID 
disordered lattice made up of Helmholtz resonators connected to a main waveg- 
uide. The disorder is introduced through the elementary cell length. An exact 
expression of the transmission coefficient of the lattice is established in the form 
of a recursive relation. Thanks to the properties of the normal distributed dis- 
order, an analytical expression of the transmission coefficient modulus and of 
the localization length is proposed. In section III, the results of this analytical 
model are compared with the results of Monte Carlo simulations based on the 
Transfer Matrix Method and a discussion is held. 

2 Propagation of acoustic waves in a weakly Gaus- 
sian disordered one-dimensional lattice 

2.1 General study 

A one-dimensional lattice is considered made up of an infinitely long cylindrical 
waveguide (with section S) connected to an array of Helmholtz resonators at 
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position Zn (Fig. [U . The Helmholtz resonators are connected to the cyhndrical 
duct through a pinpoint connection, the radius r of the throat's cross sectional 
area s of the resonators being assumed to be small compared with the wavelength 
A, of the acoustic wave {yJr/X <^ 1). 



s 
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Figure 1: Schematic representation of the Helmholtz resonator lattice made up 
of a cylindrical waveguide with section S connected to an array of Helmholtz 
resonator. 



2.1.1 Model of Helmholtz resonator 

The Helmholtz resonators are composed of a neck with section s and length i 
connected to a volume Vq (Fig. [2]). A simple model of the Helmholtz resonator 
requires the following assumptions: i) the pressure inside the volume Vb is spa- 
tially uniform, ii) the fluid in the neck moves like a solid piston. In this case, 
the air enclosed in the resonator acts as a spring for the lumped mass of air 
moving within the neck. We can furthermore consider the loss in the volume Vq 
by means of a dashpot. In these conditions, the relative change of the pressure 
p{t) — pe^^ in the volume Vb due to a small displacement x{t) of the air in the 
neck induces a restoring force F{t) 

F{t)=p{t)s^-^-^x{t), (1) 

where p is the air density and c the sound velocity |14J. The spring-dashpot 
force is considered here as linear. For a monochromatic wave (with angular 
frequency Co), the displacement x{t) =ji:e-'™ of the air is related to the acoustic 
velocity in the neck v//(f) = ve-'™ by the relation v//(f) = j(Oct:(f) and the Euler 
relation applied to the air mass m — p£eS (where £e is the effective neck length) 
submitted to the harmonic force /?/(p^e)e-'"' gives 

0)2 

jm + av + -^v = p/{pee), (2) 

where COo = cy^pp- is the eigenfrequency of the Helmholtz resonator and a 
represents the loss in the resonator. 

2.1.2 Equation propagation in the Helmholtz resonators lattice 

For a monochromatic acoustic wave with a frequency below the cut-off frequency 
of the waveguide, the acoustic pressure and velocity along the waveguide are 
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Figure 2: Elementary cell n made up of an Helmholtz resonator and a part of 
the waveguide. 



respectively written p{z,t) — p{z)e-'"" and v{z,t) — v{z)e-'"" . The amplitude p{z) 
and v(z) are related by an impedance relation. At each connection between the 
waveguide and a resonator, the wave impedance and the acoustic velocity are 
discontinuous. For Zn-i < z < Zn (Fig. [T]), the pressure and the acoustic velocity 
are denoted by p„ (z) and v„ (z) and for z„ < z < , the pressure and the acoustic 
velocity are denoted by p„+i{z) and v„+i(z). The acoustic pressure is written as 
a linear conbination of forward and backward waves (Fig. ^ as 

p„ iz) = p+e^'^-'^ + p-e--''^^--^ , (3) 

where k is the wave number defined as k = (£>/c. The continuity of the acoustic 
pressure at z = z„ yields 

Pt + Pn =Pt+\+Pn+\=P- (4) 
Moreover, the conservation of the mass flux at z = Zn leads to 

«^ = — {pt+\ - Pn+l + Pn - Pt)- (5) 

Using the Eqs. (O, ^ and ([5]), a relation between the amplitude of waves across 
the junction Zn can be written with a matrix formalism : 



Pn+l V 1 - ^ \ Pt_ 

Pn+l J \ ^ 1+W y \ p„ 



(6) 



where B = ^ and w — , ■, — — . 

We consider here the n''' cell (Fig. [21) made up of a resonator connected to the 
middle of a pipe of length t/„ and we set P^j^y, Pn ■, Pn+\ ^"^^ the amplitudes 
of forward and backward waves at the two open ends of the n''' cell. Then, the 
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relation between the pressure at z = z„ + 1/„/2 and the pressure at z = z„ — c/„/2 
can be expressed as 

where h„ = e^^'"^" . 

The ID lattice is made up of elementary cells embedded in an infinite 
waveguide. The amplitude of the incident wave on the lattice is noted the 
amplitude of the reflected wave, Pr, and the amplitude of the transmitted wave 
Pt- Using the symmetry of the lattice, the relation describing the propagation 
of a monochromatic acoustic wave through the lattice is consequently written 
as 

j =MiM2 ...M„...Mn(^ ° y (8) 

According to the Furstenberg's theorem [34j concerning the product of random 
matrix, Eq. ([5]) describes the propagation in a random media (each matrix M„ 
is different) showing localization phenomenon [55] , 

2.2 Recursive relation for the transmission coefficient 

The reflexion coefficient of the lattice Rn and transmission coefficient 7}v being 
defined by 

Pr Pt 

Rn = — and Tn = — , (9) 



the relation (IHl) is now written 



Then, if we define the matrix 5Vf as 



I I =MiM2...Mn[ ° ). (10) 



the relation 



V )=^( ? ) (12) 



and the Eq. pO| leads to m4 ~ ^ consequence, the relation 



1 

JjV-1 



( ° 1 , (13) 



allows us to express m-^ as 



m3 = w 



(1 +vv) 

hN-\ Tn-2 Tn-\ 



(14) 
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Finally, using 

/ En 
1 



MMn-iMn[ \ J, (15) 



a recursive relation between Tn, Tm^i and Tm^ can be written 



1 



Tn-1 Tn-2 ^ ^ 



Thus, the transmission coefficient Tn can be calculated from Tq = 1 and Ti = 

hi , ... rrrrn 

. In the following, results of a numerical simulation of Eq. (I16p are con- 

1 +w 

sidered as "reference results". 

2.3 Analytical expression for the transmission coefficient 

We now consider a weakly normally distributed disordered lattice with the fol- 
lowing probability density function of the cell length d„ 

1 -ii^ 

P,„W = — =.-^, (17) 



where d is the mean value of the cell length and is the variance. Regarding 
the values of d and a used in the simulation in section [3l we consider here that 
(x e] - oo; 0] ) ~ 0. If we set 

h\h2...hN 

UN^ , (18) 

■'A' 

the discrete recursive relation (|16p can then be rewritten as 

UN = [hN-\hN{l-w) + {\+w)]uN-i - [/l7V-l%]MiV-2, (19) 

with UQ — \ and u\ — \ + w. 

The value of can then be calculated using Eq. Defining h = e^^'"^ 

and using (hfj) =he^^~r- and (Z;^) = h^e^^'^''^^ , the mean value of and hpiu^ 
are 

{un) ^he-^{\-w) {liN- iUN-i) + {\ + w){uN-\) - h^e-''^''\uN-2) : (20) 

and 

{hNUN) ~ h e {\ — w){hN-\UN-\) +he [l +w){un-\) — h e ^ [un-i). 

(21) 

Finally, rewritting Eq. (pij) as follows, 

{Hnun) ^ he'^^ (un) - he'^ {1 +w){e''^~'^ - l){uN-i) , (22) 
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leads to a recursive relation concerning the mean value of as 



(UN) 



- VV) + (1 + W) - h^-^'^'"' 



(23) 

The solution of this second order recurrence equation can be obtained from the 
roots ri and r2 of the associated quadratic equation by the form below 



where 



with 



A 
B 



n 

Ci 



h{A+ ^JA^-B^ , 

= /i (a - Va2~b) , 

_ (l+w)-r2 
_ (l+w)-ri 



h{\ — w)e 



+ -il+w) 
h 



(24) 

(25) 
(26) 

(27) 
(28) 

(29) 
(30) 



Finally, the value of \Tn\ can be calculated from the expected value of upi (in 
the case of a weak Gaussian disorder) 



{\Tn\) 



1 



1 



\{un)\ iCirf + Czr; 



(31) 



The Eq. (j3ip is, in a weak disorder approximation, an analytical solution for the 
transmission coefficient modulus of a disordered ID lattice depending directly 
on the "strength" of the disorder a and on the number of cells A^. 



2.4 Analytical expression for the localization length 

In this section, the acoustic attenuation in the Helmholtz resonator is neglected, 
so that a = and w = ^^J"^^i ■ We furthermore suppose that COo < ^ which means 
that the Helmholtz resonance frequency is smaller than the Bragg frequency 
fs = Tlc/d. 

Defining rj = and 5 = ^(g-^o^i^ _ gqs. ^ and ^ become 

A = ri + (ri-;T2)5, (32) 
B = l + (2+ri2)5, (33) 
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where 
and 



Fi = cos{kd) — r| sin(A:ii) 



r2 = sin(M) +r|cos(A:ii). 
A first order expansion in 5 gives for \/A^ — B, 



Va^-b ~ ^(r? - i ) + (2r2 - 2 - - 2;rir2)5. 

Outside the stopbands of the ordered case, Fj — f < 0, and Eq. 
approximated by 



Va^-B: 



1-r? 



-rir2 + y(i-rH^) 



(34) 
(35) 

(36) 

may be 

(37) 



Lastly, this expression leads to the modulus of the roots ri and r2 as 

5. 



\ri,2\ 



2 /l T-2 



(38) 



1-rf 



Because the modulus of one of the roots is less than 1 , it can be neglected and, 
finally, the expression for the value of the transmission coefhcient outside the 
stopbands of the ordered case is 



(|rA,|) ~exp 



1 + ^L 1 



11^ 



A comparison between Eq. ([55]) and the relation 



(39) 



(40) 



giving the transmission coefficient of a disrodered lattice as a function of the 
localization length t, and the lattice length L [36], leads to the asymptotical 
expression of the adimensionnal localization length : 



t/d = lim — 

^' N^+oo\n\TN\ 



1 



11^ 



11^ 



(41) 



On one hand, the localization length is independent of the lattice length. On the 
other hand, conversely to the use of a Monte Carlo simulation, the calculation 
of the localization length ^, may be easily achieved with a low computational 
cost thank to Eq. ([4T|) . 
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3 Results and discussion 



By using the recursive relation (|16p , the propagation of an acoustic wave through 
a lattice made up of Helmholtz resonators connected to an infinite cylindrical 
waveguide with surafe area S is simulated. Each cell of the lattice has a mean 
length d. An Helmholtz resonator consists of a cylindrical volume Vq associated 
to a neck made up of a cylindrical tube of length £ and of surface area s. 

3.1 Transmission coeiRcient of a weakly Gaussian disor- 
dered ID lattice 

The Fig. ([3]) shows the well-known dispersion relation of the ordered lattice de- 
scribed above, where all the cells have the same length d. The relation dispersion 
is given by [HI 

cos{qd) = cos{kd) - ^^■?'P'?^ ^_^2^^2 sin(^^), (42) 

where q is called the Bloch wave number, Sq — s/S, Dq — d/£e and ko is de- 
fined by ko — COo/c. In the following, the lattice characteritics are chosen as 
Sq = 7.84 X 10"^ Dq = 5 and ko d ^ 0.62. 

The dispersion relation exhibits the peculiar characteristic of filters marked 
by forbidden frequencies or gaps or stopbands (marked in grey on Fig. [3]) and 
passbands in the frequency domain which result from the resonances and the 
periodic arrangements of the medium. When the relation |cos(q'ii)| < 1 is sat- 
isfied, the waves are within a passband and travel freely in the duct. On the 
contrary, when |cos(g'c/)| > 1, the waves are in a forbidden band and are spa- 
tially damped (i.e. evanescent waves). In the lattice described above, two kinds 
of stopband appear in the band structure : one is due to the resonance of the 
scatterers (Helmholtz resonators) , called resonance stopband or Helmholtz stop- 
band (marked by (a) on Fig. and the other is due to the periodicity of the 
lattice called Bragg stopband (marked by (b) on Fig. [3]). 

In Figs, m and El the moduli of the transmission coefficient of a disordered 
lattice, determined by a Monte Carlo simulation (based on the Eq. (HH])) and 
by the analytical model (from Eq. ((3T1) ). are compared. The Gaussian random 
sequence used to built the disordered lattice is characterized by a standard 
deviation a = d/50 for the Fig. [4] and a = d/20 for the Fig. [H where d is the 
mean cell length. The Monte Carlo simulation of the transmission coefficient 
is the result of 1000 realizations and is computed for a lattice made up of 200 
cells. Assuming that c/„ form an ergodic sequence, the transmission coefficient 
is also ergodic [3]. As a consequence, estimating the transmission coefficient of 
a lattice made up of elementary cells with 1000 realizations is asymptotically 
equivalent to estimating the transmission coefficient of a lattice made up of 
lOOOA^ elementary cells. The comparison of these two results shows a very good 
agreement whatever the disorder intensity. 

Firstly, the location in the frequency domain of the stopbands determined with 
the numerical and analytical methods is in good agreement with the ordered case 
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Figure 3: Dispersion relation of the ordered lattice (Eq. [42|) . The grey region 
show the stopbands of the lattice : resonance stopband (a) and Bragg stopband 
(b). 




kd 



Figure 4: Modulus of the transmission coefhcient of a weakly Gaussian disor- 
dered lattice for a = d/5Q : - analytical calculus (from Eq. (|3T|) ). □ simulation 
(from Eq. (HH)). 



result. The widths of the stopband in the disordered cases (defined for |r| close 
to zero) are increasing with the disorder intensity and, as expected, the disorder 
on the cell length clearly acts on the Bragg stopbands width [37j- Around 
the Helmholtz stopband (for 0.3 < kd < 1.1) and for a = d/5Q the maximum of 
the difference between the transmission coefficient modulus estimated by Monte 
Carlo simulation and calculated with the analytical model is 3 x lO^'*. For 
a = d/20, the maximmn of this difference is 1.6 x 10^^. 

Secondly, around the Bragg stopband, the analytical method describes very 
well the propagation of the wave in the lattice whatever the disorder intensity. 
Indeed, the comparison between both curves, for Figs. |4] and \5\ reveals only 
a weak difference for the Bragg stopband with a — d/2Q (Fig. [5|). The Fig. [6] 
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Figure 5: Modulus of the transmission coefRcient of a weakly Gaussian disor- 
dered lattice for o = d/20 : - analytical calculus (from Eq. □ simulation 
(from Eq. ([HI)). 



shows a zoom around the Bragg stopband (for 2.7 <kd < 3.6) of the transmission 
coefficient for two different disorder intensities {a — d/5Q and O = d/20). For 
the low disorder case (a — d/5Q), the difference between the analytical method 
and the simulation results for the transmission coefficient is always smaller than 
0.02. Nevertheless, with the increase of the disorder strength (for o = d/2Q), the 
comparison of the two methods shows the validity limit of the analytical model 
since the maximum of the difference between modulus of analytical transmission 
coefRcient and simulated one reaches 0.05. 
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Figure 6: (a) Modulus of the transmission coefficient of a weakly Gaussian 
disordered lattice for a = d/50 and for 2.7 < < 3.6 : - analytical calculus 
(from Eq. (I5T|) ). □ simulation (from Eq. (fTB)) V (b) Modulus of the transmission 
coefRcient of a weakly Gaussian disordered lattice for a = d/2Q and for 2.7 < 
kd < 3.6 : - analytical calculus (from Eq. pTjl ). □ simulation (from Eq. (flBl) '). 



According to Eq. PT|) . the modulus of the transmission coefRcient of a weakly 
Gaussian disorder lattice can then be easily estimated from a simple relation 
depending on the number of the cells, on the physical properties of the lattice 
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and on the standard deviation of the disorder. The computational cost is very 
low for this new method, since this calculation does not require any matrix 
product which avoids any divergence problem [37j . 

3.2 Determination of the localization length 

The analytical model presented in section [2] is used to estimate the localization 
length of a disordered lattice using Eq. (PT|) for different disorder intensities and 
for different frequencies. The lattice under consideration in this study consists 
of 2000 cells with a mean length d. The relevant adimcnsional parameter v is 
considered here for measuring the level of disorder, such that 

6a 

As a consequence, a disorder level corresponding to V = 0.5 leads to about 50 % 
of error on the lattice cell length. 

First, the analytical expression of the localization length (Eq. ([1T|) ) is compared 
to Monte Carlo simulation results (from Eq. (fTB]) ) and to the analytical results 
of Eq. (PT|) . For this, the modulus of the transmission coefficient is calculated 
using Eq. (fTB]) and (pij) and the localization length is estimated using Eq. ([^0]) . 
Figs. ([7^) and ([7)d) show this comparison for two different frequencies : (a) 
kd — 3.01, i.e., in the low edge of the Bragg stopband, (b) kd — 0.856, i.e., in 
the Hclmholtz stopband. For each case, two lattices are tested corresponding to 
N = 500 and to = 2000 cells for kd = 3.01 and to = 2000 and to = 5000 
cells for kd ~ 0.856. For each frequency, the condition — 1 < is verified. 
For kd = 0.856 the agreement between the different results may be considered 
as very good. For = 5000 and A^ = 2000, the Monte Carlo simulation and the 
calculation of the localization length with the help of the Eq. (pij) give the same 
results. The asymptotic limit (for A^ +oo) of these two methods corresponds 
to the analytical curve of the localization length ([1T|) . 

For kd = 3.01, the same remark can be made. The agreement between the 
different results is also very good. The effect of the number of lattice cells 
on the localization length estimation is clearly shown by comparing the cases 
A^ = 500 and A^ = 2000, the case A^ = 2000 joining the asymptotic limit (derived 
from Eq. (|¥T|) ) for a lower value of disorder intensity than the case oi N = 500. 

Then, Eq. (HT1) can be used to characterize the propagation in a disordered 
lattice through the localization length. Fig.[8]shows the results of such a calcula- 
tion for kd = 2.83, 2.92, 3.01, 3.03, 3.041, 3.043. AU the frequencies are in the 
low edge of the first Bragg stopband which is most influenced by the disorder 

The Fig. [S] illustrates, in a [log;log] representation, the adimensionnal lo- 
calization length as a function of the disorder strength. The influence of the 
disorder on the localization length is shown. The more disordered the lattice, 
the smaller the localization length is, whatever the frequency. Regarding the 
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Figure 7: (a) Adimensionnal localization length of a Gaussian disordered ID 
lattice vs. disorder intensity for a frequency kd ~ 0.856. (b) Adimensionnal 
localization length of a normally distributed disordered ID lattice vs. disorder 
intensity for a frequency kd = 3.01. The blue line corresponds to the analytical 
calculus (from the transmission coefRcient calculated with the help of Eq. ((3T|) ). 
the black line to the Monte Carlo simulation results (from the transmission 
coefficient estimated by Eq. (IT51) ) with 50 realizations for kd = 3.01 and 10 
realizations for kd = 0.856 and the green line to the analytical calculus of the 
localization length (from Eq. (dH)). 



frequencies, the closer they are to the center of the stopband (for kd increasing 
in the Fig. [8]), the smaller the localization length is. 

From the analytical model, it is consequently possible to evaluate the propaga- 
tion characteritics of the disordered lattice. On one hand, if the lattice length is 
smaller than the localization length, the wave is localised inside the lattice and 
the medium is considered as opaque. In the present case, for a disorder level 
V > 0.3, all the waves are localized whatever the frequencies. On the contrary, 
for V < 0.025, all the waves propagate through the lattice. 

Fig. [9] presents the adimensionnal localization length as a function of the adi- 
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Figure 8: Adimensionnal localization length of a Gaussian disordered ID lattice 
vs. disorder intensity for an adimensional frequency : (a) kd = 2.83, (b) kd = 
2.92, (c) M = 3.01, (d) kd = 3.03, (e) M = 3.041, (f) M = 3.043. 




Figure 9: Adimensionnal localization length of a Gaussian disordered ID lattice 
in function of the adimensional frequency kd for some different values of the 
disorder level : (a) 6o/<i = 0.05, (b) 6o/ii = 0.066, (c) 6o/d = 0.08, (d) 6a/d = 
0.1, (e) 6a/d = 0.13, (f) 6a/d = 0.2 and (g) 6a/d = 0.4. The frequency range is 
include in the first stopband due to Helmholtz resonance. 



mensional frequency kd for different disorder levels in the Helmholtz stopband of 
the ordered case. Contrary to above, the frequencies are chosen inside the stop- 
band in the disordered lattice case (but outside the stopband of the ordered case 
where the condition Fj — 1 < is verified). The localization lengths in this case 
are comparable with the Bragg stopband case (where the frequency are in the 
low edge of the stopband) which demonstrates that the Anderson localization 
can be different, depending on the stopband characteristics. The influence of the 
disorder level is much more important for frequencies near the Bragg stopband 
than for frequencies in the Helmholtz stopband. As a consequence, the modulus 
of the transmission coefficient is not sufficient for studying the propagation wave 
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through a disordered lattice and the analysis of the localization length can bring 
some precisions about the localization phenomenon. 

4 Conclusion 

In the present paper, a new method to study analytically the wave propagation 
in a weakly normally distributed disordered ID lattice is proposed. The trans- 
mission coefScient of the disordered lattice is calculated analytically. The results 
are in very good agreement with the Monte Carlo simulation results based on 
a recursive relation applied to the transmission coefficient. Thanks to the ana- 
lytical formulation, an expression of the localization length independant of the 
lattice length is also found. The localization length in two kinds of stopband is 
then characterized and we show that the influence of disorder on the wave prop- 
agation is more important near the Bragg stopband than near the Helmholtz 
stopband. 
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